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We have undertaken a systematic study of thermal conductivity as a function of temperature and 
magnetic field of single crystals of the compound (Lai_ z Pr z )i.2Sri.gMn207 for z(Pr) =0.2,0.4. and 
0.6. The lattice distortion due to Pr-substitution and anisotropic thermal conductivity in bilayer 
manganites are discussed on the basis of different relaxation models of local lattice distortions in 
metal and insulating states proposed by Maderda et al. The giant magnetothermal effect is scaled 
as a function of magnetization and discussed on the basis of a systematic variation of the occupation 
of the e 9 -electron orbital states due to Pr-substitution. 



I. INTRODUCTION 

Since the discovery of colossal magnetoresistance 
(CMR) in the hole-doped perovskite manganites 
Lni_2;D x MnC>3 (Ln = trivalent rare earth ; D = diva- 
lent alkaline earth Ca,Sr,Ba,Pb) , extensive studies of the 
field-induced complex physics exhibited by these CMR 
compounds have been reported^ Though the double ex- 
change(DE) mechanism 2 plays a major role in explaining 
some of the properties, a quantitative model to account 
for the insulator-metal transition and the CMR effect is 
yet to be proposed. In the DE model, the carriers are 
strongly ferromagnetically coupled, due to Hund's Rule 
coupling, to the Mn core which also results in the hop- 
ping of electrons between the Mn ions. However, it was 
pointed out that the coupling of carriers to local Jahn- 
Teller(JT) distortions and to the Mn spins must also be 
considered^. Further, several recent experiments seem 
to support the phase separation modeli. Thermal trans- 
port measurements provide additional crucial informa- 
tion on the various scattering mechanisms of thermal car- 
riers such as electrons and phonons and are expected to 
improve the understanding of some of the phenomena en- 
countered in the manganites. There have been relatively 
few reports§*£i£ on thermal transport in cubic mangan- 
ites. The thermal conductivity of cubic manganites has 
been discussed on the basis of static and dynamic lattice 
distortions by several authors. It should be noted that 
the Lni_ x D;c Mn03 cubic perovskite compounds belong 
to the Ruddlesden-Popper series generally described as 
(Ln,D) n+1 Mn „03„+i with n — oo . For n — 2, one ob- 
tains the bilayer Mn perovskite (Ln,D) 3 Mn207,in which 
two Mn02 layers are stacked with (Ln,D)2 02 layers along 
the c-axis of the structure. The reduced dimensionality 
of these compounds has been shown to have interesting 
consequences on the physical properties. Let us note that 
the optimally hole doped compound Lai.2Sri.2Mn207 



(x=0A) exhibits a CMR of more than 98% near the 
paramagnetic insulator (PI) to ferromagnetic metal (FM) 
transition temperature T c =^120K. » The thermal con- 
ductivity measurements of the double layer manganite, 
Lai.2Sri.2Mn207 were reported recently indicating an en- 
hanced giant magnetothermal effect due to its layered 
structures^. For the end member x — LaMnC>3 , which 
is a cubic compound with JT distortions, it was reported 
in that the value of k measured on a stoichiometric sin- 
gle crystal corresponds to a thermally good conductor 
, like an ideal phonon gasAi. On the other hand, the 
end member x = 1 compound , CaMn03, without JT 
active Mn sites, shows a higher thermal conduction in 
comparison with the low k data obtained for the hole 
doped compounds^. Therefore, it is desirable to examine 
a close relationship existing between the lattice distortion 
and thermal conductivity in layered manganites. 



In this paper, we report a systematic study of ther- 
mal conductivity as a function of temperature and 
magnetic field of single crystals of the compound 
(Lai_ z Pr z )i. 2 Sri. 8 Mn 2 07, for z(Pr) =0.2,0.4 and 0.6. 
Keeping the hole doping fixed at x = 0.4 , the chem- 
ical pressure effect due to substitution of Pr on La-ion 
sites enhances anisotropic lattice parameters' 2 *^ and re- 
sults in a variation of the e fl -electron character, as re- 
ported by Moritomo et ali4. The lattice distortion due 
to Pr-substitution and the anisotropic thermal conduc- 
tivity in bilayer manganites were discussed on the basis 
of different relaxation models of local lattice distortions 
in metal and insulating states proposed by Medarde et 
alii£. Next, the giant magnetothermal resistivity was 
scaled as a function of magnetization and discussed on 
the basis of a systematic variation in e g -electron orbital 
states. 
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II. EXPERIMENT 



Single crystals of (Lai-^Pr^) i. 2 Sri. 8 Mn 2 07 
(.2=0,0.2,0.4 and 0.6) were grown from sintered rods 
of the same nominal composition by the floating-zone 
method, using a mirror furnace. Crystals could be easily 
cleaved to yield shiny surfaces. X-ray Laue patterns have 
indicated the cleaved surface to be along the a6-plane of 
the structure and therefore the c-axis perpendicular to 
it. An X-ray powder pattern taken on a small part of 
the cleaved crystal did not indicate the presence of any 
additional phases. The calculated lattice parameters are 
listed in previous worksi 2 *^. The Pr-substitution at La 
ion sites resulted in a small contraction of the a-axis 
parameter and an expansion of the c-axis parameter. 
Energy-dispersive X-ray analyses of the cleaved surface 
revealed the compositions of the two samples to be close 
to those of the respective sintered rods. The dimensions 
of the crystals used in the experiments reported here 
were typically 3.4x 2.8mm 2 in the a6-plane and 1.0 mm 
along the c-axis. 

The end member compound S^M^Oy was prepared 
by a powder sintering technique under the annealing con- 
dition 1650 °C for 12hours first reported by Mizutani 
et al.ji^ and by Mitchell et al.ji*. After rapidly cool- 
ing it down to room temperatures to prevent decompo- 
sition, the polycrystalline sample was post-annealed at 
400 °C for 24hours to remove its oxygen-deficiency. The 
dimensions of the cylindrical sintered sample used in the 
thermal measurement were 4.4mm in diameter and 8mm 
along the longitudinal direction. The X-ray powder pat- 
tern of this sample showed that the sample obtained con- 
sists mainly of a single-phase with the 327 crystal struc- 
ture. The magnetization measurement also revealed that 
the present sample shows an anti-ferromagnetic transi- 
tion around T/v=~160K 17 . The thermal conductivity 
,both in the a6-plane and along the c-axis of the crys- 
tal, was measured by means of the steady-state heat-flow 
method using a G-M type cryocooler in zero field. Mag- 
netothermal conductivity measurements were carried out 
at the National Research Institute for Metals and at the 
High Field Laboratory for Superconducting Materials, 
Institute for Materials Research, Tohoku University. The 
resistivity in the a6-plane of the crystal was measured us- 
ing a conventional four-probe technique. The zero-field- 
cooled (ZFC) magnetization measurements were made 
using a SQUID magnetometer both at Iwate University 
and in Orsay. 



III. RESULTS AND DISCUSSION 

Figures Q (a) and (b) present the a6-plane normalized 
resistivity p a b and the aft-plane magnetization M a t, of our 
samples (z = 0, 0.2, 0.4 and 0.6) as a function of tempera- 
ture. The p data of the polycrystalline sample Sr 3 Mn 2 07 
are also shown in the inset of Fig.l (a). For the z — 
crystal, the value of p a b dropped sharply at ^120 K, close 




T(K) 



FIG. 1: (a)aft-plane normalized resistivity p a b and (b) aft- 
plane magnetization M a b for (Lai_ z Pr 2 ) i.2Sri.gMn207 (z = 
0, 0.2, 0.4 and 0.6) single crystals. The arrows denote the PI 
to FM transition temperature. The p data of the polycrys- 
talline sample Sr3Mn2 07 are also shown in the inset of Fig.l 
(a). 



to the insulating to metallic transition ,with decreasing 
temperature. This temperature was reduced to ^60K 
in the sample with z = 0.4. Finally, the z — 0.6 sam- 
ple showed neither an insulator to metal transition nor 
a paramagnetic to ferromagnetic transition. The tem- 
perature variation of the resistivity of the x = 1 sample 
showed a negative curvature over the whole measured 
range and the value of p rapidly increased by more than 
four orders of mag nitude, up to 2 x 10 6 ficm at 20K 
from ~100 f2cm at room temperature (RT) , which is in 
sharp contrast with the typical value of ~0.1 ficm at RT 
for Pr-substituted samplesi 2 ^. Recently , a theoretical 
band calculation 1 ^ predicts that the bilayer perovskite 
SraM^Cv is an antiferromagnetic insulator with an in- 
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FIG. 2: (a)afr-plane and (b)c-axis thermal conductivity of 
(Lai_ z Pr 2 ) i.2Sri. 8 Mn 2 7 (z=0, 0.2 ,0.4 and 0.6) single crys- 
tals. The arrows denote the PI to FM transition temperature. 
For comparison, the temperature variation of poly crystalline 
Sr3Mn207 is also presented in Fig. 2(a). 



direct gap of 0.45 eV , which is qualitatively in good 
agreement with our magnetization and resistivity data. 
Moreover, the electronic structure of this compound is 
expected to be similar to that of the cubic perovskite 
end member CaMnC>3. 

The afe-plane and the c-axis thermal conductivity of 
(Lai_ z Pr z ) i. 2 Sri. 8 Mn 2 7 (z=0, 0.2 ,0.4 and 0.6) single 
crystals are shown in Fig|2]as a function of temperature, 
up to 300K. The arrows denote the PI to FM transition 
temperature. For comparison, the temperature variation 
of the end member polycrystalline sample Sr3Mn207 is 
also presented in Fig. 2 (a). We note, first that a low ther- 
mal conductivity was observed over a wide of tempera- 
ture range for all samples ,in spite of the high-quality of 
our single-crystals. Second, the temperature dependence 



of k in the afe-plane exhibited anomalies associated with 
the PI to FM transition, while no such anomalies were 
observed for k along the c-axis direction near Tc for any 
of the samples. Above Tc, the value of K a b showed a pos- 
itive curvature but tended towards a constant value near 
room temperatures with increasing temperatures. The 
electron thermal conductivity ,for the z = sample ,is 
estimated to be at most ~1 % of the total value, even 
in the metallic state from the resistivity data based on 
the Wiedemann-Franz (W-F) law and reproduce, neither 
the whole anomaly in n a b, nor the giant thermal effect 
around Ic^fi. Thus, heat carriers are taken as phonons 
for all measured samples. Here, the phonon mean free 
path l p h in the insulating state of the x = 0.4 sample is 
estimated from the specific heat data C p h, and longitu- 
dinal sound velocity v p h , using the kinetic expression for 
the phonon gas n — 1/3 • C p h,v p hl P h- If C p / l (150K)=1.78 
J/cm 3 Kii2> and « o &(150K)=11.3 mW/cm (present data) 
for the (Lai_ a Pr z ) i. 2 Sri. 8 Mn 2 07 single crystal (z = 0.4) 
and u J) / l =4.3xl0 5 cm/sSi in the afr-plane for the La 
i.2Sri.gMn 2 07 bilayer single crystal , we then obtain 
lph,ab =4.43 A at 150K which is comparable with the Mn- 
O-Mn bond length (~4 A). In a smilar way, the value of 
l p h,ab reaches 4.07 A at 300K. In our estimation, the value 
of v a b is taken as a constant since the relative variation 
of the longitudinal velocity is very small from ultrasonic 
measurements on the bilayer manganites. 

Next, the end member compound Sr3Mn 2 07 showed a 
different behavior from the afe-plane data in the com- 
pound (Lai_^Pr 2 ) i. 2 Sri. 8 Mn 2 07. The value of n of 
the x = 1 sample remained almost a constant down 
to ~ 200K, then showed a slight increase upon crossing 
Tn =~ 160 K and finally dropped rapidly with decreas- 
ing temperature. The temperature dependence of k of 
the x = 1 sample reminds us of a typical phonon conduc- 
tion in insulating materials, except for the magnitude of 
measured values. The extremely high resistivity of the 
x = l(n = 2) sample ( ~100 Ocm at RT) in sharp con- 
trast with the value of p a b (~ O.lfkm at RT) for the 
Pr-substituted x = 0.4(rt = 2) samples gives rise to a 
negligible contribution of small polarons ,which proba- 
bly results in a typical phonon conduction. In compari- 
son with the thermal data of the polycrystalline sample 
of Sr3Mn 2 C>7 , low thermal values of single crystal sam- 
ples are anomalous indicating a strong phonon damping 
which correlates with dynamical lattice distortions cou- 
pled with thermally hopping of localized carriers at Mn- 
sites (as discussed later). 

We discuss the relationship between local lattice dis- 
tortions and thermal conductivity in bilayer manganites. 
Figures |3] (a) and (b) show the anisotropic lattice stric- 
tion and volume striction for the z=0, 0.2 and 0.4 sam- 
ples, respectively. The anomaly in anisotropic lattice 
striction reverses sign near z=0.2 but the volume of all 
samples rapidly shrinks upon crossing Tc- These find- 
ings have been discussed in ref.21 on the basis of differ- 
ent relaxation models of local lattice distortions in metal 
and insulating states proposed by Maderda et alii^. The 
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FIG. 3: (a)anisotropic lattice striction and (b)volume stric- 
tion,as a function of temperature, up to 300K, for (Lai_ z Pr z ) 
i.2Sri. 8 Mn 2 07 (z = 0, 0.2, 0.4 and 0.6) single crystals. The 
arrows denote the PI to FM transition temperature. 



authors tried to explain the different behaviors of av- 
erage and local lattice distortions observed in neutron 
powder diffraction experiments on La 2 _ 2 a;Sri + 2 ;i; Mn207 
(0.32<x<0.4). Based on their explanation, the observed 
volume shrinkage near Tq seems to arise from a strong 
suppression in the local lattice disorder manifested by 
Mn 3+ -0 6 and Mn 4+ -0 6 octahedra with the JT distorted 
and undistorted sites when the system is going towards 
the metallic state from the insulating state. The MnO 
bond-length distribution becomes narrower below Tc be- 
cause of a screening effect due to itinerant carriers. Ac- 
cordingly , the phonon conduction is limited by the lo- 
cal lattice disorder accompanying a polaron hopping in 
the insulating state. However, in the metallic region , 
it is expected that the homogenous structure of the lo- 
cal lattice associated with carrier delocalization, causes 



an upturn in the afr-plane thermal conductivity, below 
Tc ■ It should be noted that this conclusion is reached 
from the lattice striction and thermal data of the same 
samples, combined with the preceding model. Further- 
more, in comparison with the anisotropic Debye- Waller 
(D-W) factor of Sr 3 Mn 2 07 , the anisotropic D-W data 
of La2-2xSri_|_2a;Mn207 (0.32<x<0.4) show that, even in 
the metallic state of the latter, the inhomogeneities in the 
local lattice structure still remain. This finding is prob- 
ably related, not only with a small upturn in k, but also 
to the low values of K. On the other hand, the out-of- 
plane phonon conduction across the MnO double-layers 
is strongly disturbed from stacking fault-like scattering 
, of the layered structures. Phonon scattering due to 
stacking faults dominates thermal conduction along the 
c-axis, masks some contribution of the local lattice dis- 
tortion associated with the insulator to metal transition 
and results in to the anisotropy of the thermal conduc- 
tivity. 

We discuss next, the normalized magnetothermal con- 
ductivity (NMTC) [n ab (H)/k & (0)] of the Pr-substituted 
crystals at z=0.2 , z=0A and z=0.6 (Figgjl. The field 
was in the ab-plane and parallel to the direction of the 
heat current. Our data reveal several interesting features. 
First, a giant magnetothermal conductivity was observed 
in association with the occurrence of the CMR effect and 
the NMTC reached about 30%, up to 8T, for all sam- 
ples in spite of a lower Tc- Here, it should be noted 
that k c (H)/ k c (0) is almost negligible in contrast with a 
rapid enhancement of K a b (H)/K a {, (0), which is consistent 
with the absence of anomalies observed in k c (T)^. Sec- 
ond, accompanying the field-induced first-order PI to FM 
transition of the insulating sample, z=0.6 , a clear hys- 
teresis in the magnetothermal conductivity was observed. 
The remnant NMTC at 4.2K reached nearly the same 
value of ~10% ,both in first and second scans, where, 
after the first scan, the sample was warmed up to high 
temperatures in order to demagnetize it. The value of 
NMTC of the z — 0.6 sample exhibited a step-like behav- 
ior at 4.2 K, while it showed a smooth rise with increasing 
temperature. As previously pointed out, electronic con- 
tribution from the W-F law estimation using the CMR 
data reproduces neither an upturn in K a b nor a giant 
thermal effect. Magnetostriction data around Tc showed 
a volume shrinkage upon application of a magnetic field, 
up to 8T, as well as spontaneous striction (~0.1%) upon 
crossing Tp ?2? . The negative magnetovolume effect indi- 
cates a suppression in local lattice disorder ;the reduction 
in phonon scattering due to structural disorder leads to 
a giant magnetothermal effect in the same way as the 
upturn in k below Tc- 

The normalized magnetothermal resistivity 
[Kab(M) I « o &(0)] -1 of our samples is plotted as a 
function of magnetization at selected temperatures in 
Fig. 5 (a,b,c), where the insets in Figs. 5 (a) , 5 (b) 
and 5 (c)represent M a y-H. curves at the corresponding 
temperatures for the samples with 2=0, 0.2 and 0.6 
respectively. It is easily found that the value of the 
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FIG. 4: normalized magnetothermal conductivity (NMTC) 
Kb (tt)/K ab (0)] of (Lai_ 2 Pr z ) i. 2 Sri. 8 Mn 2 7 (z = 0, 0.2, 
0.4 and 0.6) single crystals.(a)2=0.2;(b) 2=0.4 and (c) 2=0.6. 
The field was applied in the ab-plane, parallel to the direction 
of the heat current. 



magnetothermal resistivity is scaled with a parabolic 
function of 1 — A(M/M S ) 2 ,as pointed out for cubic 
manganites , by Cohn et al£. Here, M s is taken as the 
value of 3.6fiB per Mn-sites. The coefficient A depends 
on the value of z and shows a monotonical decrease from 
0.4 for z=0 down to 0.28 for z=0.6, as shown in Fig5.(d). 
This finding suggests a suppression of the lattice-spin 
coupling in thermal CMR due to Pr-substitution and 
can be explained as follows. The Pr-substitution had 
a strong effect on the magnetic anisotropics . The 
low-temperature magnetic anisotropy M a j,/M c in a field 
of lOOmT reduced from ~10 at 2=0 to nearly 1 at z=0.2 
and 0.4. For z=0.6, the easy axis of magnetization was 
along the c-axis. On the other hand, the lattice parame- 
ters vary with Pr-substitution in such a manner that the 



FIG. 5: Normalized magnetothermal resistivity 

[K a b(M)/ KabiO)]^ 1 of (Lai_ 2 Pr z ) i. 2 Sri.8Mn 2 7 {z = 0, 
0.2, 0.4 and 0.6) single crystals as a function of magnetiza- 
tion, at selected temperature. (a)2=0;(b) 2=0.2 and (c) 2=0.4 
and 0.6. The insets in Figs.5(a,b,c)represent M a b-H curves 
at the corresponding temperatures for the samples with 2=0 
, 0.2 and 0.6 respectively, (d) z dependence of coefficient A 
and lattice parameter c/a 



c-axis expands but the a(b) axis shrinks (see Fig. 5 (d)) 
which results in a change of e g -electron character from 
planar 3d x 2_ y 2 to linear 3g? 3z 2_,,2. Accordingly, such a 
variation of the e ff -electron character accompanies that 
of magnetic anisotropy, so that Pr-substitution probably 
suppresses indirectly the lattice-spin coupling in thermal 
CMR through a variation of spin alignment, strongly 
coupled with orbital states. In other words, a higher 
magnetothermal conduction favors the occurrence of 
the easy axis of magnetization lying in the a6-plane. 
It is an interesting feature that a change of e s -electron 
character ,due to Pr-substitution, appears in giant 
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,z=0.2 crystal was observed. If we accept the above ex- 
planation, this finding indicates that phonon scattering 
due to local lattice disorder ( bond-length fluctuation) 
remains, even in the metallic state, which is consistent 
with the anisotropic D-W data of La2-2a;Sri + 2 a ;Mn207 
(0.32<x<0.4) in the metallic state. Recently, Maezono et 
al^i pointed out that in manganite systems, the dynamic 
JT effect is more enhanced in the metallic region than in 
the insulating region because of strong electron correla- 
tions. Magnetothermal measurements ,at low tempera- 
tures , probe the lattice disorder not only in the insu- 
lating region but also in the metallic region as the D-W 
factor does. Furthermore, it is desirable to clarify the 
crucial role played by the dynamic local lattice distor- 
tion on phonon conduction ,not only in the insulating 
phase , but also in the metallic phase through magne- 
tothermal measurements on high quality single crystal of 
the extremely doped bilayer compounds. 

In summary, a systematic study of thermal conductiv- 
ity as a function of temperature and magnetic field of sin- 
gle crystals of the compound (Lai_aPr z ) 12 Sr 1 . 8 Mn207 
for z(Pr) =0.2,0.4. and 0.6 was performed. The a6-plane 
phonon conduction of the Pr-substituted x = 0.4 samples 
is disturbed by the dynamical lattice distortions associ- 
ated with polaron hopping on Mn sites. On the other 
hands, the typical phonon behavior observed for the end 
member x=l sample is probably attributed to a reduc- 
tion of the inhomogeneities in the local lattice distortions 
strongly coupled with small polarons. The giant magne- 
tothermal effect has been explained taking into account a 
systematic variation of the e 9 -electron orbital states due 
to Pr-substitution. 



magnetothermal effect involving heat carrying phonons. 

Finally, we shall stress that at lower temperatures, far 
from Tc, a finite increase of K a h ( ~14% at 8T) for the 



Acknowledgments 

The authors would like to thank Dr. H.Ogasawara for 
his technical support. This work was supported by a 
Grant-in- Aid for Scientific Research from the Ministry of 
Education, Science and Culture, Japan. 



* Electronic address: matsukawa@iwate-u.ac.jp 

1 Colossal Magnetoresistance, Charge Ordering and Related 
properties of Manganese Oxides (Ed. C. N. R. Rao and B. 
Raveau, World Scientific, Singapore, 1988); Colossal mag- 
netoresistive oxides (Ed Y. Tokura, Gordon and Breach 
Science. Publisher, 2000). 

2 C.Zener, Phys. Rev. 82, 403 (1951); P.G.de Gennes ibid. 
118, 141 (1960). 

3 A. J. Millis, P.B. Littlewood, and B.I. Shraiman, Phys. 
Rev. Lett. 74, 5144 (1995); A.J. Millis, B.I. Shraiman,and 
R.Mueller, ibid. 77, 175 (1996). 

4 For a recent review, see E. Dagotto,T. Hotta and A. Moreo, 
Phys. Rep.344, 1 (2001). 

5 D.W. Visser, A. P. Ramirez, and M.A. Subramanian, Phys. 
Rev. Lett. 78, 3947 (1997). 



6 J.L. Cohn, J. J. Neumeier, C. P. Popoviciu, K. J. McClel- 
lan, and Th. Leventouri, Phys. Rev. B 56, R8495 (1997). 

7 S. Uhlenbruck, B. Buchner, R. Gross, A. Freimuth, A. 
Maria de Leon Guevara and A. Revcolevschi, Phys. Rev. 
B 57, R5571 (1998). 

8 J. Hejtmanek, Z. Jirak, M. Marysko,S. Krupicka, C. Mar- 
tin , A. Maignan, M. Hervieu and B. Raveau, Phys. Rev. B 
60, 14057 (1999). 

9 Y. Moritomo, A. Asamitsu, H. Kuwahara and Y. Tokura, 
Nature, 380, 141 (1996). 

M.Matsukawa , H. Ogasawara, R. Sato , M. Yoshizawa , 
R. Suryanarayanan , G. Dhalenne , A. Revcolevschi and 
K. Itoh, Phys. Rev. B62, 5327 (2000). 

1 J.S.Zhou and J.B. Goodenough, Phys. Rev. B64, 024421-1 
(2001). 



7 



H.Ogasawara, M.Matsukawa, S.Hatakeyama, 

M.Yoshizawa, M.Apostu, R.Suryanarayanan, G.Dhalenne, 
A.Revcolevschi, K.Itoh and N.Kobayashi, J. Phys. Soc. 
Jpn. 69 1274 (2000). 

M.Apostu, R.Suryanarayanan, A.Revcolevschi, 

H.Ogasawara, M.Matsukawa, M.Yoshizawa and 

N.Kobayashi, Phys. Rev. B64 (2001) (012407). 

Y. Moritomo, Y. Maruyama, T. Akimoto and A. Naka- 

mura : Phys. Rev. B56 (1997) R7057. 

M.Medarde, J.F. Mitchell, J.E. Millburn, S. Short and J. 

D. Jorgensen: Phys. Rev. Lett. 83 (1999) 1223. 

N.Mizutani,A.Kitazawa,O.Nobuyuki and M.Kato, 

J.Chem.Soc.(Japan)Ind.Ed.73(1970)1097. 

J. F. Mitchell, J.E. Millburn, M.Medarde, D.N. Argyriou 



and J. D. Jorgensen, J.Appl.Phys.85(1999)4352. 

18 H.Meskine,Z.S.Popovic and S.Satpathy, Phys.Rev.B65 
(2002) 094402. 

19 A.Vasilev,T.Voloshok, M.Apostu, R.Suryanarayanan, 
A.Revcolevschi ,JETP Lett. 73 (2001) 630. 

20 C.N.Brosseau,M.Poirier, R.Suryanarayanan, G.Dhalenne, 
A.Revcolevschi, cond- mat/0012380 

21 M.Matsukawa, H.Ogasawara,T. Sasaki, M.Yoshizawa, 
M.Apostu, R.Suryanarayanan, A.Revcolevschi, K.Itoh 
and N.Kobayashi, J. Phys. Soc. Jpn. 71 (2002)1475 . 

22 T.Kimura,Y.Tomioka,A.Asamitsu and Y.Tokura, 
Phys.Rev.Lett. 81 (1998) 5920. 

23 R.Maezono and N.Nagaosa cond-mat/0202301 



